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Abstract. We rediscuss tiie temperature of O dwarfs based on new non-LTE lin e blan lteted atmosphere models 
including stellar winds computed with the CMFGEN code of Hillier & Miller ( 1998| ). Compared to the latest 
calibration of Vacca et al. (199C), the inclusion of line blanketing leads to lower effective temperatures, typically 
by ~ 4000 to 1500 K for 03 to 09.5 dwarf stars. The dependence of the Teg-scale on stellar and model parameters 
- such as mass loss, microturbulence, and metallicity - is explored, and model predictions are compared to optical 
observations of O stars. Even for an SMC metallicity we find a non-negligible effect of line blanketing on the T^s- 
scale. The temperature reduction implies downward revisions of luminosities by ~ 0.1 dex and Lyman continuum 
fluxes Qo by approximately 40% for dwarfs of a given spectral type. 
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As a significant fraction of the flux of O stars is emitted 
in the inaccessible Lyman continuum (A < 912 A) reli- 
able direct determinations of their effective temperatures 
are not possible. Indirect methods, primarily based on at- 
mospheric modeling, are therefore employed (e.g. Bohm- 
Vitense 1981 , Crowther 1997 ). Given the need for a de- 
tailed treatment of non-LTE effects a nd th e presence of 
stellar winds (Kudritzki & Hummer 1990|) , a complete 
modeling of such atmospheres including also the effects of 
numerous metal-lines ( "line blanketing" ) remains a com- 
plex task (cf. Schaerer fc Schmutz |1994 Hillier & Miller 



1998| , Pauldrach et al. ^OOlD 



For these reasons, most published spectral analysis 
have so far been based on simple non-LTE models. For 
example, the most recent calibration of stellar parame- 
ters of O and early B type stars of Vacca et al. (1996, 
hereafter VGS96) is based only on results from plane par- 
allel, pure Hydrogen and Helium (H-He) non-LTE models. 
Their derived temperature scale for O stars is found to be 
significantly hotter than most earlier calibrations (see ref- 
erences in VGS96). Such differences lead to non- negligible 
changes in the fundamental parameters of O stars — e.g. 
luminosities, Lyman continuum fluxes etc. — when esti- 
mated from spectral types. Accurate calibrations are cru- 
cial for various astrophysical topics, such as comparisons 
with stellar evolution models, determinations of the initial 
mass function and cluster ages, studies of H ii regions, and 
others. 



Indications for a decrease of Toff due to line blanket- 
ing effects have been found since the flrst non-LTE + 
wind modeling attempts by Abbott & Hummer ( |1985| , 
and subsequent investigations based on the same "wind 
blanketed" models), the improved models of Schaerer 



& Schmutz ( |1994D and Schmutz ( |1998D , and the fully- 
blanketed plane parallel non-LTE models of Hubeny et 
al. (1998). Similar indications are obtained by FuUerton 
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et al. (200C) from recent modeling of FUSE spectra with 
t he co de of Pauldrach et al. ( |2001| ) and by Crowther et al. 
( ^00 1|) . 

The effective temperature scale of O stars is revised 
here ba sed o n the recent CMFGEN code of Hillier & 
Miller (|l998|) , which treats the problem of a non-LTE 
line blanketed atmosphere with a stellar wind in a direct 
way, thereby avoiding possible shortcomings due to opac- 
ity sampling techniques employed by Schaerer & Schmutz 
( |1994| ), Schmutz ( [l998[ ), and Pauldrach et al. ( |200l| ). First 
results on the dwarf sequence are presented here. A more 
detailed account including all luminosity classes will be 
presented in a subsequent publication. 

In Sect. 2 we describe our method and the calculated 
models. The results, their dependence on model/stellar 
parameters, and first comparison with observations are 
presented in Sect. 3. Implications of the revised T^s scale 
and remaining uncertainties are discussed in Sect. 4. 



2. Model ingredients 

We have constructed spherically expanding non-LTE 
line-blanketed model atmospheres using the CMFGEN 
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comovme -frame code of Hillier & Miller ( |199S| ). This code 
solves the equations of statistical equilibrium, radiative 
transfer, and radiative equilibrium, and allows for a direct 
treatment of line blanketing through the use of a super- 
level approach. The following ions are included in our cal- 
culations: H, He i-ii, C ii-iv, N ii-v, O ii-vi, Si ii-iv, 
S IV- VI, and Fe iii-vii, whose ~ 2000 levels are described 
by ^ 700 super-levels, corresponding to a total of ~ 20000 
bound-bound transitions. 

For simplicity a constant Doppler profile (thermal 
width corresponding to the mass of Helium and T — 20000 
K plus a microturbulcnt velocity of wturb — 20 km/s) is as- 
sumed for all lines in the statistical equilibrium and radia- 
tive transfer computation. To examine if a constant ther- 
mal width and the use of the large microturbulcnt velocity 
does not artificially enhance the photospheric blanketing, 
we have made test calculations with the correct depth and 
ion dependent thermal width and i^turb = 0.1 km/s. No 
significant changes in atmospheric structure, level popu- 
lations, and the emergent spectrum were found. This is 
explained in part by the high density of lines in the UV 
part of the spectrum, which implies an average spacing 
between lines which is smaller than the typical Doppler 
width. The opacity in the wing of a line is therefore mostly 
dominated by the core opacity of the neighbouring line, 
and the exact intrinsic line profile is of little importance. 
With our standard choice, ~ 80000 frequency points are 
necessary to correctly sample all lines. 

The input atmospheric structure, connecting smoothly 
the spherically extended hydrostatic layers with the wind 
(parametrised by the usual /3-law), is calculated as in 
Schaerer & de Kot er ([1997D with the ISA-WIND code 
of de Koter et al. (1996) As the approximate tempera- 
ture structure in ISA-WIND differs from the final radia- 
tive equilibrium temperature structure, the atmosphere 
structure in the quasi-hydrostatic part may be inconsis- 
tent with the final gas pressure gradient. However, for the 
issues discussed here the differences are small (correspond- 
ing to a change of ^ 0.1 dex in logg). In any case, the lines 
considered here are formed in the transition region whose 
structure/dynamics remain largely parametrised. The for- 
mal solution of the radiative transfer equation yielding the 
detailed emergent spectrum allows for incoherent electron 
scattering and includes standard Stark broadening tables 
for H, He i, and He ii lines. Our standard calculations 
assume Wtmb = 5 km/s. 

We have computed a grid of models representative 
of O dwarfs in the temperature range between ^ 30000 
and 50000 K. The model parameters are taken from the 
CoStar models A2-E2 of Schaerer & de Koter ( |1997D , with 
an additional model Y2 at {Tes,\ogg) - (31500,4.0) and 
the remaining parameters |^ taken from stellar tracks of 
Meynet et al. (1994). For each parameter set a line blan- 
keted model with solar metallicity and a pure H-He model 
was computed. 



3. Results 

3.1. Blanketing effect on the temperature scale 

The optical He i A4471 and He ii A4542 classification 
lines are used to assign spectral types to our models. 
Fig. shows the effective temperature as a function of 
logW' = logVF(4471) -logW(4542) and the correspond- 
ing spectral type according to Mathys ( 19881 ). The pure 
H-He models (open circles) follow closely the Tcff-scale for 
dwarfs of VGS96, which is based on a compilation of stel- 
lar parameters determined using pure H-He plane parallel 
non-LTE model atmospheres. The comparison shows that 
if we neglect line blanketing our dwarf model grid would 
yield nearly the same absolute Tcff-scale as the pure H- 
He plane parallel models adopted for the spectral analysis 
included in the compilation of VGS96. 

The line blanketed model sequence (Fig. |^, filled sym- 
bols) shows a systematic shift to earlier spectral types for 
a given temperature, or equivalently a shift to lower Tes 
for line blanketed models at a given spectral type. The dif- 
ference ranges from ~ 1500 K at spectral type 09.5 to 
~4000 K at spectral type 03 (cf. Fig. |l], solid line in lower 
panel). The difference with the VGS96 scale is shown as 
the dotted line. Our line blanketed scale smoothly joins 
earlier cahbrations at 09.7V (see VGS96, Fig. 1). 

As a spectral type corresponds to a given ionisation 
state of Helium in the line formation region, blanketed 
models must be more ionised than unblanketed models. 
The introduction of line blanketing leads to three main ef- 
fects illustrated in Fig. || for the c ase of model C2 (cf. Figs. 
13 and 14 of Schaerer & Schmutz 1994). Qualitatively the 
same trends are obtained for all models. 

1 ) Blanketing leads to the backscattering of photons to- 
wards the inner atmosphere which forces the local tem- 
perature to rise so that flux conservation is fulfilled 
(backwarming effect; see upper panel). 

2) At the same time the radiation field becomes more 
diffuse, as quantified by the dilution factor W — 1 — 
\F/ J shown in the middle panel, causing an in crease 
of the mean intensi ty (cf . Abbott & Hummer 1985 , 
Schaerer & Schmutz 1994). 

3) In the outer part of the atmosphere (logr^oss ^ ~2 
in Fig. H) the ionisation is essentially controlled by 
the EUV flux, which is quite strongly reduced due to 
the blocking by numerous metal lines shown in Fig. 
0. Here this effect dominates over 2), in contrast with 



^ M = 16.83 M0, logTeff = 4.498, log(L/LQ) = 4.552, R = 
6.358 Rq, logM = -7.204 Mq/jv, and = 2500 km/s. 



the finding of Schaerer & Schmutz (1994), leading to 
a lower ionisation. 

Effects 1) and 2) lead to a higher ionisation in the 
formation region of the classification lines. This results 
predominantly in an increase of W(4542) at Toff ^ 38000 
K and a decrease of W(4471) at higher Toff(cf. Fig. |). 

Given the stronger mass loss and the corresponding 
increase of the wind density, one expects even larger tem- 
perature differences between non-blanketed and line blan- 
keted models for giant and supergiant luminosity classes 




Fig. 1. Upper panel: Effective temperature of O dwarfs 
as a function of the spectral subtype (lower scale). The 
correspondance between log W (upper scale) and spectral 
type is given by Mathys (1988). For values logVF' > 1.0 
we assign a spectral type of 10 Filled circles show our line 
blanketed models, open circles pure H-Hc models. The 
VGS96 relation (dotted line) is well reproduced by our 
pure H-He models. Lower panel: T^g shift between H-He 
and line blanketed models (solid line) and between VGS96 
scale and our line blanketed models (dotted line). Note the 
decrease of T^g due to line blanketing. 



Fig. 2. Comparison of atmosphere structures of model 
C2 (Tcff=41.8 kK, log 5 =4.0). Solid line is for the line- 
blanketed model and dashed line for the pure H-He 
model. Upper panel: Temperature structure. Middle panel: 
Dilution factor W = 1 — jF/J where F is the flux and 
J the frequency averaged mean intensity (cf. Schaerer & 
Schmutz 1994 ). Lower panel: Populations of the ground 



levels of Helium and of the lower and upper levels of the 
transitions He i A4471 and He ii A4542 . Given are the rel- 
ative number population n^ with respect to total H popu- 
lation ntot(H). 



(cf. Abb ott & Hummer |l985i Schmutz |l998| , Crowther et 



al. POOID 



3.2. Dependence on model and stellar parameters 

How strongly do our results depend on poorly known pa- 
rameters such as the velocity law in the photosphere-wind 
transition zone, Wturbj and variations of gravity and M ex- 
pected within the dwarf class ? Do our calculations still 
miss opacity sources ? 



As pointed out by Schaerer & Schmutz (1994) changes 
in He line profiles due to modifications of the velocity 
law v{r) in the photosphere-wind transition zone can lead 
to similar equivalent widths variations as line blanketing. 
Test calculations for models A2 and C2 varying the slope 
(3 from 0.8 (our standard value) to 1.5 show that both 
H-He and line blanketed models exhibit a similar shift 
in log(Vl^'). The obtained relative T^g difference between 
H-He and blanketed models remains thus identical. The 
blanketed models with /3 = 1.5 have log(M^') lowered by 
~ 0.1-0.2 dex. However, as Ha profile fits for O dwarfs 
are generally quite compatible with f3 ^ 0.8 (e.g. Puis et 



al. 1996 ), we do not expect drastic changes of the absolute 
scale from this effect. 

An increase of the microturbulent velocity Wturb from 
5 to 20 km/s in blanketed models increases the strength 
of He I A4471 (cf. Smith & Howarth Villamariz & 

Herrero 200C), and leads to a shift of + 0.05 to 0.1 dex 
in log(W^') (i.e. towards later types) for models with TcB 
Si 42000 K. For hotter stars the difference is negligible. 

The effect of line blanketing is strengthened further in 



denser winds (cf. Abbott & Hummer 1985, Schmutz 1998). 
Models C2 and D2 with an increased mass loss rate by a 
factor of 2 show a shift of log(W^') between ~ - 0.05 and 
- 0.1 dex. 

Test calculations for model C2 including also Nickel 
(Ni IV- vi) show unchanged He lines. Other models in- 
cluding also Ar, Ne, and Ca confirm that Fe blanketing 
dominates. 

While microturbulence and mass loss affect (though in 
opposite ways) the exact Teff-scale, their exact importance 
will have to be studied in future comparisons. 
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Fig. 3. UV spectrum of model C2 with line blanketing 
(solid line) and pure H-He model (dashed line). Note the 
reduction of the EUV flux below ^ 500 A due to the in- 
clusion of metals. 

3.3. Comparison with observations 

As a first comparison of our models with observations 
we show in Fig. ^ the predicted and observed equiva- 
lent widths of He i and He ii classification lines and other 
strong He lines frequently used in spectral analysis. The 
observational data is taken from Mathys ( |1984 |1989D and 
Conti & Alschuler (1971). The observational scatter is 
real, as the typical measurement errors are 5-7 %. The 
general trend is that the He i A4471 and He i A4388 equiva- 
lent widths are well represented by the models, while He ii 
A4542 seems to be overestimated by 20% for spectral 
types earlier than 07. He ii A4200 behaves as He ii A4542 . 
The other equivalent widths remain essentially unchanged 
by all other parameter variations discussed above (Sect. 
3.2). A value of (3 <; 1.5, a stronger increase of M, or an 
unrealisticly large reduction of log g would be necessary to 
reduce the predicted equivalent widths of the Stark broad- 
ened He II lines. 

Strictly speaking, if we were simply to reduce VF(4542) 
by ~ 20% while keeping M^(4471) constant for early spec- 
tral types, this would result in a change of ~ - 0.08 dex 
in log(W^') thus reducing the shift in the Tcff-scale be- 
tween line blanketed and pure H-Hc models from ~ 4000 
K to 3000 K in the high temperature part. Future tailored 
spectral analysis should allow to assess more precisely the 
achievable fit accuracy and th e pr ecise importance of the 
parameters discussed in Sect. 3.2 on the stellar parame- 
ters. 
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Fig. 4. Comparison between observed (filled squares: lu- 
minosity class V; open squares: other luminosity classes) 
and calculated equivalent widths of He i A4471, He i 
A4388, He ii A4542 , and He ii A4200 (in A). Line blanketed 
models are indicated by full circles, pure H-He models by 
open circles. See discussion in text. 



3.4. Comparison with previous analysis 

As discussed in Sect. 1, few earlier studies have addressed 
the effect of line blanketing in O stars. Essentially all inves- 
tigations concur with a reduction of Toff when blanketing 
is included. 

Abbott & Hummer ( 1985 ) have constructed a core- 



halo model where backscattered radiation due to multiple 
line scattering in the wind modifies the plane parallel pho- 
tosphere. Their so-called "wind blanketed" models yield a 
decrease of Toff by ~ 10 % for 04 types (similar to our re- 
sults), ~ - 2000 K for an 09.5 supergiant, but e ssentially 
no shift for 09.5 dwarfs (Bohannan et al. 1990 , Voels et 
al. 1989). The latter finding is likely due to lack of photo- 
spheric blanketing (inherent to their method) and modest 
wind blanketing due to the comparatively low mass loss 
rates of 09.5 dwarfs. 

An improved Monte-Carlo opacity sampling method of 
a unified photosphere-wind model was used by Schaerer 



& Schmutz ( |1994| ), Schaerer & de Koter ( |1997D , and sub- 
sequently applied to a larger parameter space by Schmutz 
(1998). For mass loss rates comparable to the values 



adopted here (typica l for dwarfs with low mass loss) the 
models of Schmutz ( |l998|) indicate differences from ^ - 

600 K at 08 to 2000 K at 04, which is half the 

shift deduced from Fig. |^ and roughly the difference ob- 
tained with Z = 1/8 Z© (see Sect. 0). This indicates that 
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their method underestimates hne blanketing compared to 
CMFGEN. 

Using plane parallel line blanketed non-LTE models 
based on opacity distribution functions Hubeny, Heap & 



Lanz (199S) found that a pure H-He model with T^ff 
37500 K and log 3=4.0 is necessary to reproduce the H 
and He lines of a line blanketed model with Tcs = 35000 
K and same gravity. As can be seen from Fig. ^ our results 
are in excellent agreement with their result. 

LTE line blocking has bee n inc luded in plane paral- 
lel models by Herrero et al. ( 2000 ) primarily to resolve 
discrepancies between He i singlet and triplet lines. For 
stars with Tes ^ 40000 K this leads to a strengthening 
of He I A4471, opposite to the effect found in all above 
studies including ours. This results must be due to an in- 
complete treatment of the various effects of line blanketing 
(cf. above), and appears to be unphysical. This discrep- 
ancy with line blanketed models has also been noted by 
the authors. 



4. Implications and concluding remarks 

The importance of line blanketing obviously depends on 
metallicity Z . Therefore one may wonder at which Z the 
stellar parameters will again correspond to the results ob- 
tained with pure H-He (metal-free) atmosphere models, 
i.e. close to the VGS96 scale. Test calculations for models 
A2 and D2 with a metallicity close to the SMC value (1/8 
Z0) show still a reduction of Toff compared to pure H-He 
models : A Tcff is ^ 60 % that found at solar metallicity. 

As the bolometric correction is essentially unchanged 
by line blanketing, and the Aly versus spectral type 
(Sp) calibration independent to first order from model- 
ing, we can use the BC-Tcff relation of VGS96 to derive 
luminosities through log(L/Lo) = 2.736 log Tcs{Sp) — 
AMv{Sp) —9.164:. This relation shows that the predicted 
reduction of T^s by ^ 0.04 dex implies a downward revi- 
sion of L by ^ 0.1 dex for dwarfs of a given spectral type. 

Since line blanketing is mostly efficient in the EUV, 
the ionising spectrum below 912 A is modified. The to- 
tal number of Lyman continuum photons Qo predicted 
by our models is in good agreement with the calcula- 
tions of Schaerer & de Koter (1997). The change of Qo 
due to the shift in the Tcff-Sp calibration, taking into 
account the change of both the radius and the ionising 
flux per unit surface area qo, is given by A(log(5o) = 
-1.264A(logTeff) + A(loggo(rcff)), wh ere th e latter term 
is dominant (sec Schaerer & de Koter 1997). For a given 



spectral type between 04V and 09V this amounts typi- 
cally to a reduction of Qo by ~ 40 %. 

While the results presented here provide a clear im- 
provement over earlier calibrations, and a general reduc- 
tion of Tcff due to line blanketing is unavoidable, we wish 
to caution that the absolute Tog scale may still be subject 
to revisions for the following reasons. First, tailored multi- 
wavelength analysis of individual objects are required to 
test the present models in mor e dep th for O stars, as re - 
cently started by Bouret et al. ( |2GGlD , HiUicr et al. ( |200lD , 



and Crowther et al. ( ^001 ). Second, the effect of X-rays 
on the overall ionisation balance and in particular on the 
Helium lines remains to be studied. Indeed for late O and 
B stars, depending on the relative X-ray to photospheric 
flux at energies close to the relevant ionisation poten- 
tials and the wind density. X-ray emission (likely due to 
shocks) is expected to increase the ionisation of most ions 
(MacFarlane et al. 1994). Nonetheless, first test calcula- 
tions with CMFGEN seem to indicate that photospheric 
lines are not affected by X-rays generated in the wind. 
Finally, we note that comparisons of photoionisation mod- 
els calculated using fluxes from recent atmosphere models 
(including CMFGEN and Pauldrach et al. |200l| models) 
with ISO observations of H 11 regions possibly reveal a flux 



deficiency at energies ^ 34.8-40.9 eV (Morisset et al. 2001 



but cf. Giveon et al. 2001 ). The importance of the latter 
two findings — possibly related to each other — on the 
lines used here as Teg indicators remains to be studied. 

As UV and optical classification lines of O stars de- 
pend in fact on several parameters {Tcs, gravity, mass loss 
rate, met allicit y, rotation; e.g. Abbott & Hummer |1985, 
Walborn et al 



Schmutz 1998 



1995), spectral type and 



luminosity class calibrations must ultimately account for 
this multi dimensionality. Some of these issues will be ad- 
dressed in subsequent publications. 
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